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\ Abstract: In gravity-mediated SUSY breaking models with non- universal gaugino masses, 

' lowering the 5'f7(3) gaugino mass IM3I leads to a reduction in the squark and gluino masses. 

' Lower third generation squark masses, in turn, diminish the effect of a large top quark 

Yukawa coupling in the running of the higgs mass parameter rn?fj^, leading to a reduction 
in the magnitude of the superpotential // parameter (relative to Mi and M2). A low 
r~| , 1^1 parameter gives rise to mixed higgsino dark matter (MHDM), which can efficiently 

annihilate in the early universe to give a dark matter relic density in accord with WMAP 
^ \ measurements. We explore the phenomenology of the low IM3I scenario, and find for the 

5— ( ' case of MHDM increased rates for direct and indirect detection of neutralino dark matter 



relative to the mSUGRA model. The sparticle mass spectrum is characterized by relatively 
light gluinos, frequently with rrig <C m-g. If scalar masses are large, then gluinos can be very 
light, with g — > Zig loop decays dominating the gluino branching fraction. Top squarks 
can be much lighter than sbottom and first/second generation squarks. The presence of 
low mass higgsino-like charginos and neutralinos is expected at the CERN LHC. The small 
— mass gap should give rise to a visible opposite-sign/same flavor dilepton mass 
edge. At a TeV scale linear e^e~ collider, the region of MHDM will mean that the entire 
spectrum of charginos and neutralinos are amongst the lightest sparticles, and are most 
likely to be produced at observable rates, allowing for a complete reconstruction of the 
gaugino- higgsino sector. 
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1. Introduction 



Recently, a variety of astrophysical measurements by the WMAP[Q] and other collabora- 
tions have determined the density of cold dark matter (CDM) in the universe to be 

ncDAih"^ = 0.113 ±0.009. (1.1) 

The additional determination of a non-zero dark energy component to the universe suggests 
that we live in a ACDM universe, with f^A/i^ ~ 0.35. While the nature of dark energy 
remains a mystery, there are a number of well-motivated particle physics candidates for 
the CDM, and collider and DM search experiments may serve to distinguish between the 
various possibilities in the near future. 

One of the especially intriguing features of i?-parity conserving supersymmetric models 
is that they provide a natural candidate for cold dark matter (CDM) in the universe |2|, 
|3|. The lightest neutralino Zi in gravity-mediated SUSY breaking models (SUGRA) is 
especially appealing as a DM candidate in that it can be produced thermally in the early 
universe, and the calculable relic abundance turns out to be in the right neighborhood to 
match the measurements of the density of CDM in the universe. 

Many analyses of neutralino CDM have been performed Q within the context of the 
paradigm minimal supergravity model|p (mSUGRA), which is completely specified by the 
parameter set 

mo, mi/25 Aq, tan/3 and sign{^). (1-2) 

The mSUGRA model assumes that the minimal supersymmetric Standard Model (MSSM) 
is valid between the mass scales Q = Mqut and Q = M^eak- A common value mo 
(mi/2) ((^0)) is assumed for all scalar mass (gaugino mass) ((trilinear soft SUSY breaking)) 
parameters at Q = Mgut, and tan/? is the ratio of vacuum expectation values of the two 
Higgs fields that give masses to the up and down type fermions. The magnitude of the 
superpotential Higgs mass term /i, but not its sign, is fixed so as to reproduce the observed 
Z boson mass. The values of couplings and other model parameters renormalized at the 
weak scale can be computed via renormalization group (RG) evolution from Q = Mqux 
to Q = M„cak- From these weak scale parameters, sparticle masses and mixings may be 
computed, and the associated relic density of neutralinos as well as scattering cross sections 
and decay rates can be determined. 

In most of the allowed mSUGRA parameter space, the relic density ^^_^h'^ turns out 
to be considerably larger than the WMAP value. Consistency with WMAP thus implies 
that neutralinos should be able to annihilate very efficiently in the early universe. In 
the mSUGRA model, the annihilation rate is enhanced in just the following regions of 
parameter space, where the sparticle masses and/or the neutralino composition assume 
special forms. 

• The bulk region occurs at low values of mo and m^/2[@) In this region, neutralino 
annihilation is enhanced by t-channel exchange of relatively light sleptons. The bulk 
region, featured prominently in many early analyses of the relic density, has been 
squeezed from below by the LEP2 bound on the chargino mass mfr, > 103.5 GeV 
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and the measured value of the branching fraction B(b — > 57), and from above by the 
tight bound from WMAP. 

• The stau co-annihilation region occurs at low mo for almost any mi/2 value where 

~ The staus, being charged, can annihilate rapidly so that fiZi co- 

annihilation processes that maintain Zi in thermal equilibrium with fi, serve to 
reduce the relic density of neutralinos Q . 

• The hyperbolic branch/focus point (HB/FP) region at large mo ~ several TeV, 
where |//| becomes small, and neutralinos efficiently annihilate via their higgsino 
components 1^]. This is the case of mixed higgsino dark matter (MHDM). 

• The ^-annihilation funnel occurs at large tan f3 values when 2m^^ ~ m^i (or mn) 
and neutralinos can efficiently annihilate through the relatively broad A and H Higgs 
resonances IQ]. 

In addition, a less prominent light Higgs h annihilation corridor occurs at low mi/2|lC] and 



a top squark co-annihilation region occurs at particular Aq values when m^^ ~ m^^\ll\. 

Many analyses have also been performed for gravity-mediated SUSY breaking models 
with non-universal soft SUSY breaking (SSB) terms. Non-universality of SSB scalar masses 
can, 1. pull one or more scalar masses to low values so that "bulk" annihilation via f-channel 



exchange of light scalars can occur ||T^, 13|, 2. bring in new near degeneracies of various 



sparticles with the Zi so that new co-annihilation regions open up [H, ll,!!, 3. bring the 
value of rriA into accord with 2m^ so that Higgs resonance annihilation can occur |l3|| , 
or 4. pull the value of down so that higgsino annihilation can occur [p^, [T^ . 

If non-universal gaugino masses are allowed, then qualitatively new possibilities arise 



that are not realized in the mSUGRA model|18, |2^, 21 1. One case, that of mixed 
wino dark matter (MWDM), has been addressed in a previous paper ||22|. In this case, 
as the weak scale value of SU{2) gaugino mass M2(weak) is lowered from its mSUGRA 
value, keeping the hypercharge gaugino mass Afi(weak) fixed, the wino component of Zi 
continuously increases until it becomes dominant when M2(weak) < Mi (weak) (assuming 
is large). The ZiWi^2W coupling becomes large when Zi becomes wino-like, resulting 
in enhanced ZiZi W~^W~ annihilations. Moreover, co-annihilations with the lightest 
chargino and with the next-to-lightest neutralino help to further suppress the thermal relic 
abundance of the lightest SUSY particles (LSPs). Indeed, if the wino component of the 
neutralino is too large, this annihilation rate is very big and the neutralino relic density 
falls well below the WMAP value. 

A qualitatively different case arises in supersymmetric models if the SSB gaugino 
masses Mi and M2 are of opposite sign. In this case, as \Mi\ and IM2I approach one 
another, there is little bino-wino mixing, and the Z^ maintains a nearly pure bino-like or 
wino-like identity. The WMAP relic density can nonetheless be achieved for Mi ~ —M2 via 
bino-wino co- annihilation (BWCA) of the bino-like lightest neutralino. The resulting DM 



and collider phenomenology was investigated in Ref. |23]. The MWDM and BWCA DM 



scenarios were also investigated recently in Ref. |24], where these scenarios were collectively 
dubbed "the well-tempered neutralino". 
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In this paper, we investigate a scenario where, as previously noted by Belanger et al.\(25\ 



and also by Nezri and Mambrini[26|, a diminution of the GUT scale value of the SU{3) 
gaugino mass M3 relative to Mi and M2 leads to a sparticle spectrum with lower gluino 
and squark masses (the latter are lowered through RG effects due to a reduced M3). We are 
motivated to consider this because by adjusting M3 to the right value(s) provides another 
way to obtain MHDM. To understand this, we begin by noting that the RG equation for 
the soft SUSY breaking Higgs squared mass fnjj^ reads (in a standard notation [p7|]) 

where Xt = mQ^+m'^ +m\^+A^ and S = m?^^ — m|^^+Tr(mg — m| — 2m^ + m|)+m|;). 
Here, ft is the top quark Yukawa coupling and t = logQ^. The ffXt term drives "rn^fj^ 
to negative values due to the large top quark Yukawa coupling in the celebrated radiative 
electroweak symmetry breaking (REWSB) mechanism [p8| . In the case where M3 << Mi ~ 
M2 at the GUT scale, the squark squared mass terms and (and hence Xt) are suppressed 
at lower scales; as a consequence, rri^j^ is not driven to such large negative values as in the 
universal gaugino mass case. Thus, if IM3I is chosen small (but not so small that m-^^ is 
no longer driven negative), we still obtain REWSB, but with a smaller weak scale value of 
—mjf^} There is also a corresponding effect on the RG flow of it^h^, but this is typically 
smaller because <C ft except for very large values of tan /3. The weak scale value of 
(at tree-level) is then obtained from the weak scale parameters of the Higgs sector via the 
EWSB relation, 



~ (tan2/3-l) - ■ ^'-^^ 



We see that if ^ and moderate to large values of tan/3, fj? ~ —m?^^. Thus the 

smaller In-z^^J value expected in the low IM3I case results in a smaller |//| parameter, and 
a correspondingly larger higgsino component of the lightest neutralino Zi. 

This situation is illustrated in Fig. |l|, where we plot the evolution of fn^^, ^Hu^ 
and the gaugino mass parameters versus the renormalization scale Q from Q = Mqut to 
Q = M„eak for the mSUGRA model with niQ = mi/2 = 300 GeV, Aq = 0, tan /? = 10, > 
and nit = 175 GeV (solid curves), and for the case of M3 = 0.5mi/2 (dashed curves). The 
electroweak gaugino mass parameters evolve identically at the one loop level, and the tiny 
difference seen is a two loop effect. The gluino mass parameter, on the other hand, starts 
off at a smaller value and evolves to a correspondingly smaller value at the weak scale. 
Turning to the mass parameters in the Higgs sector, we see that, as expected, m'jj^ runs 
to a less negative value in the case of the low M3 model than in the case of the mSUGRA 
model with universality of the GUT scale gaugino mass parameters. The evolution of m'jj^ 
is very similar in the two cases because the bottom Yukawa coupling is very small. 

This reduction of is illustrated as well in Fig. |2|, where we show the behaviour of /j, as 
a function of the ratio \M3\/mi/2, for several values of mo, again with Ao = 0, tan/3 = 10, 
H > and rrit = 175 GeV. The curves end where REWSB is no longer viable because rnjj 



""^ Indeed, for given values of other parameters, the constraint of REWSB imposes a lower bound on jA/a 
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m„ =300GeV, m^^ =300GeV, tanp=10, Ag =0, n >0, m, =175GeV 
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Figure 1: Evolution of the soft SUSY brealcing Higgs mass parameters sign(TO|^ 



and 



sign {m'jj^)^\m'j^J as a function of scale Q in the mSUGRA model (solid) for mo = 300 GeV, 
mi/2 = 300 GeV, Aq — 0, tan/3 = 10, /i > and mt — 175 GeV. The same running mass 
parameters are shown for LM3DM for the same parameters as in the mSUGRA case except taking 
M3 = 0.5mi/2 at Mgut (dashes). Also shown is the corresponding evolution of gaugino mass 
parameters Mi, M2 and M3, for the mSUGRA case (solid) and the M3 — Q.bmi/2 case (dashes). 



does not evolve to negative values. We see that for relatively low values of mg (smaller 
than a few times mi 12), low values of ^ are achieved for \M^\/mii2 < 1, while for very large 
values of (that, for the chosen value of mi/2 may have been forbidden for the mSUGRA 
case), REWSB with low values of |^| becomes possible but only for \M3\/mi/2 > 1- These 
low regions are just generalizations of the well-known HB/FP regions of the mSUGRA 
model. The location of the "generalized" HP/FB region in the rriQ — m\j2 plane of the 
extended model depends on the value of rj, = M^/mii2'- it lies to the left (right) of the 
corresponding region in the mSUGRA model if |r3| < 1 (|r3| > 1). In the study presented 
here we focus on the first possibility, as it leads to lighter coloured sparticles that may 
well be accessible at the CERN Large Hadron Collider (LHC) scheduled to commence 
operations next year. At the same time, the large higgsino component expected in this 
low M3 dark matter (LM3DM) scenario should lead to larger detection rates relative to 
mSUGRA in direct and indirect searches for neutralino dark matter. 

Many previous studies have examined the neutralino relic density in models with gaug- 
ino mass non-universality, along with prospects for direct and indirect detection of DM 
neutralinos. Griest and Roszkowski|29] first pointed out that a wide range of relic density 
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Figure 2: The values of /i dictated by REWSB as a function of |M3|/mi/2i for Ao = 0, tan f3 — 10, 
/Lt > and mt = 175 GeV, at various values of mo 



values could be obtained by abandoning gaugino mass universality by allowing departures 
from M1/M2 — 0.5. Anomaly-mediated SUSY breaking models, where the gaugino masses 
are proportional to the /3-functions of the corresponding low energy gauge groups have 
Ml : M2 : M3 ~ 3 : 1 : — 10. As a result, the Zi is almost a pure wino which annihilates 
very efficiently, resulting in too low a thermal relic neutralino density: to account for the 
observed dark matter density, Moroi and Randall |]30[| invoked the decay of heavy moduli to 
wino-like neutralinos in the early history of the universe. The dark matter relic density and 
detection rates in models with non-minimal SU (5) gauge kinetic function, and also in O-II 
string models were studied by Corsetti and Nath[^]. Birkedal-Hanson and Nelson showed 
that a GUT scale ratio M1/M2 ~ 1.5 would bring the relic density into accord with the 
measured CDM density via MWDM, and also presented direct detection rates [ 32 1. Berlin, 
Nezri and Orloff studied the variation of relic density and the enhancements in direct and 
indirect DM detection rates as gaugino mass ratios are varied] 33]. Bottino et al. performed 
scans over independent weak scale parameters to show variation in indirect DM detection 
rates, and noted that neutralinos as low as 6 GeV are allowed|34]. Mambrini and Muhoz, 
and also Cerdeno and Muiioz, examined direct and indirect detection rates for models 
with non- universal scalar and gaugino masses] 35]. Auto et a/.] 14] proposed non-universal 
gaugino masses to reconcile the predicted relic density in models with Yukawa coupling 
unification with the WMAP result. Masiero, Profumo and Ullio exhibit the relic density 
and direct and indirect detection rates in split supersymmetry where Mi, M2 and are 
taken as independent weak scale parameters with ultra-heavy squarks and sleptons[^]. Fi- 
nally, as mentioned above, the variation of the relic density due to the change of M3 - the 
subject of this paper - was first studied by Belanger et al. who showed that large swaths 
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of the mo — mi/2 plane are consistent with the WMAP value when the SU{3) gaugino 
mass M3 becomes small |25|; this topic was subsequently also studied by Mambrini and 
Nezrill]. 

It has been shown that the various non-universal scenarios each lead to distinctive 
phenomenologies, and can be distinguished from mSUGRA and from one another via 
their implications for accelerator experiments, and simultaneously, for direct and indirect 
searches for DM. The purpose of this paper is to study WMAP viable SUSY models with 
a non- universal GUT scale gaugino mass hierarchy IM3I <C Mi ~ M2 - these models have 
received relatively little attention in the literature - and to explore their phenomenology. 
In regions of parameter space that yield the observed relic density of MHDM, we examine 
prospects for its direct and indirect detection, and also outline the impact on prospects 
for direct detection of sparticles at the Fermilab Tevatron, the CERN LHC and at the 
future international linear e~^e~ collider (ILC). For expediency, we adopt an mSUGRA- 
like model with universal GUT scale SSB parameters, but with the SU{3) gaugino mass 
as one additional parameter; i.e. we take Mi = M2 = rni^2 > at Q = Mqut, while 
allowing M^{Mgut) to remain as a free parameter with either sign. By dialing IM3I to 
low enough values (for the range of m-o that we consider) any point in the remainder of 
the parameter space can be WMAP allowed. The parameter space naturally divides into 
regions with bino dark matter (EDM), or with MHDM. Once the WMAP constraint is 
fulfilled, then in the MHDM case one finds generally enhanced rates for direct and indirect 
DM detection. As far as colliders go, a mass spectrum with rriq ~ is predicted in the 
scalar sector. In the gaugino sector, a much reduced mass gap of rrig — m^^^ is expected as 
compared to mSUGRA. This means in part that lighter gluinos can be allowed despite the 
constraints from LEP2, and that the Fermilab Tevatron may explore a substantial portion 
of the LM3DM parameter space via gluino pair production. We find that in the portion 
of the parameter space where mg/M2 is most suppressed, tjiq is necessarily large, and the 
radiative decays g — > gZi constitute the dominant decay modes of the gluino. In this case, 
gluino pair production may lead to dijet events at hadron colliders. At the CERN 

LHC, an enhanced reach is found in tuq vs. m\j2 parameter space relative to the mSUGRA 
model due to the reduced squark and gluino masses. At a linear e^e~ collider, a much 
lighter spectrum of squarks and gluinos is expected. In the case of MHDM, the low [i 
parameter implies that the entire spectrum of charginos and neutralinos is rather light, 
and may be accessible to ILC searches for new sparticles. 

The remainder of this paper is organized as follows. In Sec. we outline the parameter 
space of the LM3DM scenario, and show that any point in parameter space may be WMAP 
allowed if a suitably low value of IM3I is adopted. We also illustrate characteristic features 
of the sparticle mass spectrum expected in this scenario. In Sec. ^, we discuss expectations 
for direct and indirect detection of neutralino DM in the LM3DM scenario, and show that 
generally enhanced detection rates are expected if MHDM occurs. In Sec. ^ we give an 
overview of some of the main features of the LM3DM scenario which give rise to distinct 
signatures at the Fermilab Tevatron, the CERN LHC and a = 0.5 — 1 TeV ILC. In Sec. 
Ill, we present a summary of our results and some conclusions. 

Note added: As this paper was being finalized, the WMAP collaboration released their 
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three year data along with implications for cosmology 1 37]. Their new central value of 
^^CDM^^ (and of course, the quoted error) is slightly lower than in However, the 

central value obtained by combining the WMAP data with other data as in Table 6 of Ref. 
1 37] is almost unchanged from ( |1.1| ). In either case, our analysis is hardly affected. 

2. Parameter space and mass spectrum in the LM3DM scenario 

As discussed in the previous section, the LM3DM scenario is completely specified by the 
parameter set: 

mo, Aq, tan/3, and sign{fj,), (2-1) 

(together with rrit which we fix to be 175 GeV), where we assume Mi = M2 = ^ 
at Q = Mgut, and where M3 can assume either sign. The assumed equality of Mi and 
M2 can be relaxed somewhat and our conclusions suffer little qualitative change so long 
as Ml ~ M2. To calculate the sparticle mass spectrum, we adopt Isajet 7.73[^], which 
allows for the input of non-universal scalar and gaugino masses in gravity mediated SUSY 
breaking models where electroweak symmetry is broken radiatively. The relic density is 
evaluated via the IsaReD program]39], which is part of the Isatools package. IsaReD 



evaluates all 2 — > 2 tree level neutralino annihilation and co-annihilation processes and 
implements relativistic thermal averaging in the relic density calculation. 

We begin our discussion with an illustration to show that any point in mSUGRA model 
parameter space that is WMAP disallowed owing to too large a relic density Q^^h"^, can 
become WMAP allowed by adjusting {M^l until a suitably small \^\ value that yields a 
relic density in accord with WMAP is attained. Assuming that mo/m-1/2 is not too large, 
this is achieved by lowering {M^l from its mSUGRA value. As an example, we plot in Fig. 
|3|a the neutralino relic density ^^_h'^ vs. M3 for mo = mi/2 = 300 GeV, while Aq = 0, 
tan/3 = 10 and > 0. The value M3 = 300 GeV puts us in the mSUGRA model, and here 
we see ^^Jt^ = 1-1) and so the parameter space point would be excluded. As ilfs is lowered 
from its mSUGRA value, gluino and squark masses also decrease, resulting in a lower weak 
scale value of ]'7Z^^] and hence \^\ as discussed in Sec. |l|. At M3 = 150 GeV, the value 
of ]/u] has diminished sufficiently that the Zi is no longer bino-like, but is instead a mixed 

higgsino-bino state. This is illustrated in frame 6), where we plot Rfj = \l v^i^"^ + v^^"^ (in 
the notation of Ref. [^]), which gives an indication of the higgsino components of the Zi. 
As M3 is decreased even further, the relic density reaches a minimum around M3 ~ 110 
GeV, and then increases slightly before reaching the LEP2 limit where m^^ becomes less 
than 103.5 GeV. This slight increase occurs because m^^ drops below Mz and then My/-, 
so that .^1.2^1 ZZ, W~^W~ processes, the major LSP annihilation modes in the early 
universe become kinematically suppressed. Since there is no reason to favour the positive 
sign of /Vf3, we show the behavior of the relic density and for negative M3 as well, and 
note that these are nearly symmetrical under M3 —M3. 

The effect on the sparticle mass spectrum of lowering the magnitude of the SU{3) 
gaugino mass M^ is shown in Fig. ^, where we plot the sparticle mass spectrum versus the 
ratio r3 = M^/mii2 for the same parameters as in Fig. |^. At r3 = 1, we see the usual 
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mQ=300 GeV, my^=300 GeV, tanp=10, Aq=0, ^i>0, m =175 GeV 
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Figure 3: a) The neutralino relic density ilcDAih^, and b) higgsino component i?^ of the Ughtest 
neutrahno as a function of A/3 for niQ — 300 GeV, mi/2 — 300 GeV, Aq — 0, tan/3 = 10, /i > and 
TO* = 175 GeV. 

hierarchy of sparticle masses as obtained in the mSUGRA model. As M3 is lowered, the 
gluino masses reduces sharply from nig = 727 GeV in the mSUGRA case to rrig ~ 400 
GeV for = 0.5 where ^^.^h"^ = 0.11. The reduction of M3 reduces the renormalization 
suffered by the squark mass parameters, and causes the squark masses to correspondingly 
drop from the vicinity of 700 GeV in mSUGRA to ~ 500 GeV for = 0.5. On the other 
hand, slepton mass parameters, whose renormalization does not depend on SUSY QCD 
effects at the one loop level, are hardly affected by the change of M3. Thus, in the LM3DM 
scenario, the mass scale of squarks and sleptons is more nearly equal, and less hierarchical, 
than the case of mSUGRA at low rriQ. In fact, for the case shown in Fig. ^, the top squark 
ti has dropped to a lower mass than the various sleptons in the case of = 0.5. The 
value of the parameter is shown by the black dotted curve, and we see that this drops 
sharply as IM3I decreases. The drop in increases the higgsino component of the lighter 
charginos and neutralinos and, once they cross over to becoming higgsino-like, their masses 
decrease with decreasing M3 as well. 

We also show in Fig. ^ the sparticle masses for negative values of r^. The slepton and 
first/second generation squark masses are nearly symmetrical about M3 = because the 
leading (one loop) contributions to their renormalization group evolution are quadratic in 
the gaugino masses. However, the top squark and various chargino and neutralino masses 
are not symmetric. For the stops, this occurs because the RG evolution of the A parameters 
(At in this case) depends linearly on the gaugino masses. The asymmetric evolution of the 
A-parameters also affects m'jj^ and m'jj^ the same way; as a result, and hence chargino 
and neutralino masses are also not symmetric under M3 — > —M^. We see, in fact, when 
M3 < 0, the WMAP measured value of ^IcDMh"^ is attained at a value of = —0.6, in 
contrast to = 0.5 for positive values of M3. 
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LM3DM: irig =300GeV, m^,^ =300GeV, fanp =10, Ag =0, n >0, m, =175GeV 
r> 900 




Figure 4: A plot of various sparticle and Higgs boson masses and the /i parameter vs. 
for Too = 300 GeV, TOi/2 = 300 GeV, = 0, tan/3 = 10 and ^ > 0. 

Various sparticle masses are shown in Table Q corresponding to the parameters shown 
in Fig. ^. We show the spectrum for the mSUGRA case, together with that for the 
LM3DM1 case (with rs = 0.5) and the LM3DM2 case, with rs = -0.6. In the several 
rows (below the masses) we show the relic density, BF{b 57), Aa^, and the neutralino- 
proton scattering cross section a{pZi) for these scenarios. Finally, in the last row, we show 
that component of the "up-type" neutral higgsino, i.e. that couples to the T3 = +1/2 
quark-squark system, in the neutralino LSP: we will return to this, as well as the LM3DM 
scenario in the last column of Table |l| when we discuss Tevatron signals for the gluino in 
the LM3DM framework. ^ 

In Fig. ^, we show contours of in the uiq vs. mi/2 plane for Aq = 0, tan /3 = 10 and 
// > 0, where at each point in the plane IM3I has been reduced until the value of the relic 

^We are aware that the value of nth, especially in the LM3DM1 and LM3DM2 scenarios, is well below 
the limit from searches for the SM Higgs boson which should be applicable in these cases because ttia is 
large. Since the value of depends on our choice of tan/3 as well of Aq (whose precise values do not 
qualitatively affect the features that we discuss here), we will continue to use these scenarios as simple 
illustrations of the model. 
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Table 1: Masses and parameters in GeV units for mSUGRA and three LM3DM scenarios. In each 
case, ^0 = 0, tan/3 = 10 and mt = 175 GeV. 



density is found to be ^^^h^ = 0.11. Frame a) shows these contours for negative rs, while 
frame b) shows contours for r3 > 0. The red shaded region on the left hand side is excluded 

because the fi becomes the LSP (in contradiction to search hmits for stable charged or 
colored relics), while the blue region at low is excluded by the LEP2 chargino search 

limits, which require iriy^^ > 103.5 GeV. Unlike in the mSUGRA model, the LEP-excluded 
blue region is not fiat because the chargino is gaugino-like for small values of mo but 
higgsino-like as mo becomes large. Of course, in the stau co-annihilation region, which is 
tight against the boundary of the red region, the relic density is already in accord with the 
WMAP value even in the mSUGRA case. In the lower-left bulk region, IM3I need only be 
reduced to values of ~ —0.7 (frame a)) or ~ 0.6 (frame b)). However, for large niQ 
and low mi/2, values of [r^l as low as ~ 0.3 are needed to reach accord with WMAP. A 
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LM3DM: -M^ < m^,^ , tan^=10, A„ =0, [i >0, m, =175 GeV LM3DM: < m^,^ , tan^=10, A„ =0, [i >0, m, =175 GeV 
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Figure 5: Contours of a) < in the mo vs. mi/2 plane for mi/2 up to 1.5 TcV, with tan f3 = 10, 
= 0, /i > 0. Each point in the plane has dialed to such a value that ^^i^^ ~ O-H- The red 
region on the left is excluded either because fi becomes the LSP or electroweak symmetry is not 
correctly broken, while the blue region is excluded by the LEP lower limit > 103.5 GeV. In 

frame &), we plot contours of ra > for the same parameter choices as in frame a), although we 
extend the range of mi/2 to 3 TeV. 

third case study in this region, labelled LM3DM3, is shown in the last column of Table |l|.^ 
In most of parameter space, values of |r3| ~ 0.4 — 0.6 are sufficient to match the measured 
relic density. 

The planes of Fig. ^ naturally divide into two distinct regions. The left-hand side of 
each plot with |r3| ~ 0.5 — 0.6 is labelled the bino DM region (BDM), since here the Zi is 
bino-like, while the large mo side of the plane with |r3| ~ 0.5 — 0.6 is labelled as MHDM, 
since here the Zi is mixed higgsino-bino state. The bino-wino-higgsino content of Zi for 
an mi/2 = 1005 GeV slice out of Fig. |5| is shown in Fig. ^. 

In the BDM regions of Fig. ^, the values of M3 and uiq are low enough that 
approaches 2m^^. In this case, the thermal distribution of neutralinos convoluted with the 
ZiZi bb cross section allows for an enhanced annihilation rate via the s-channel A, H- 
pole. On the very low mo edge of parameter space, fi co-annihilation also contributes to a 
reduction in the neutralino relic density, and a wide range of rs is possible. As we move to 
larger niQ values in Fig. the value of niA becomes much larger than 2m^^, and A- funnel 
annihilation is no longer efficient enough to reduce the relic density. Various sparticle and 
Higgs masses are shown as a function of rriQ in Fig. ^, for the same parameters as in Fig. 
|6[ In this case, for large uiq, must be reduced to lower values of < 0.5, and the 
Zi becomes MHDM. Then ZiZi WW, ZZ, and also W\Z\ and (to a smaller extent) 

^In this case, the decays Z4 qqg and the decays W2 ~~* iJ-dg, not included in Isajet, are kinematically 
accessible. We expect, however, that this will not significantly afltect our analysis because these sparticles 
will dominantly decay via their two-body modes. 



- 11 - 



mj,2=1005 GeV, M3<m,/y tanP=10, \=0, \i>0, m =175 GeV 
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Figure 6: The bino, wino and higgsino content of the Zi versus tuq for a mi/2 — 1005 GeV shce 
out of the plane in Fig. showing that Zi is EDM for small mo, and MHDM for large mg. 

Z1Z2 co-annihilation act to suppress the relic density. This is shown explicitly in Fig. 
m, where we plot the thermally averaged neutralino annihilation cross sections integrated 
over temperature from freeze-out to the present time, versus mo for the same parameter 
choices as in Fig. ^. We see that the co-annihilation processes become significant only for 
niQ ~ 400 GeV where r^, and correspondingly also have become sufficiently small. 

Finally, we mention that in the upper-left of frame of Fig. ^6) (which has been extended 
to — 3 TeV to facilitate the discussion of the LHC reach in Sec. [4.2| ), the value 

drops below 0.5 near the contours at mo ^ 0.8 TeV and ~ 3 TeV. In this part of 

the plane, a small additional reduction in M3 is needed to match flcDAih'^ = 0.11 because 
the enhancement to ZiZi annihilation via amplitudes with s-channel H, A exchanges no 
longer obtains since 2m^ > rriH, ruA- Thus, in this upper-left region, the Zi is again 
MHDM. 




3. Direct and indirect detection of neutralino CDM 



In this section, we explore the prospects for direct and indirect detection of neutralino dark 



matter within the LM3DM framework |4C]. We adopt the DarkSUSY code [41], interfaced 
to Isajet, for the computation of the various indirect detection rates, and resort to the 



Adiabatically Contracted N03 Halo model|42| for the dark matter distribution in the Milky 
Way, which tends to give higher detection rates, especially for gamma ray and anti-particle 
detection than other halo profiles. In this respect, our projections may be regarded as 
optimistic.^ We evaluate the following neutralino DM detection rates: 



Direct neutralino detection via underground cryogenic detectors [ 46 1. Here, we com^ 



pute the spin independent neutralino-proton scattering cross section, and compare 
it to expected sensitivities [|3 for Stage 2 detectors (CDMS2||l, Edelweiss2[|9|, 



*For a comparison of the implications of different halo model choices for indirect DM detection rates, 
see e.g. Refs. ||, 0, ^, 0. 



- 12 - 



n7,/2 =1005GeV, < m^,^ , tan^=10, A„ =0, n >0, m, =175GeV 
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Figure 7: Sparticle masses vs. mo for an mi/2 = 1005 GeV slice out of the plane in Fig. ||6). 



CRESST2|g, ZEPLIN2l5|) and for Stage 3, SuperCDMS, along with ton-size de- 
tectors (XENON H, GERDAH, ZEPLIN40 and WARP|5|). We adopt the 
projected (mass-dependent) sensitivities of CDMS2 and 1-ton XENON detectors as 
the experimental benchmark for direct DM detection at stage 2 and stage 3 detectors. 

Indirect detection of neutralinos via neutralino annihilation to neutrinos in the core of 



the Sun[|56[. Here, we present rates for detection of — > conversions at Antares[|57| 
or IceCube|^^. The reference experimental sensitivity we use is that of IceCube, with 
a muon energy threshold of 50 GeV, corresponding to a flux of about 10 muons per 
km^ per year. 

• Indirect detection of neutralinos via neutrahno annihilation in the galactic center 
leading to gamma rays[^], as searched for by EGRET [Q, and in the future by 
GLAST|61|. We evaluate the integrated continuum 7 ray flux above a = 1 GeV 
threshold, and assume a GLAST sensitivity of 1.0x10" 
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• Indirect detection of neutralinos via neutralino annihilations in the galactic halo 
leading to cosmic antiparticles, including positrons]^ (HEAT||6^, Pamela []64|| and 
AMS-02|6|), antiprotons|6| (BESS|3l, Pamela, AMS-02) and anti-deuterons {D) 
(BESS|6^, AMS-02, GAPS|^). For positrons and antiprotons we evaluate the av- 
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m,,2 =1005GeV, < iriy^ , fanp=JO, =0, [i >0, m, =175GeV 
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Figure 8: The thermally averaged neutralino annihilation cross sections times relative velocity 
integrated from a; = to xf versus mo for a mi/2 = 1005 GeV slice out of the plane in Fig. 



eraged differential antiparticle flux in a projected energy bin centered at a kinetic 
energy of 20 GeV, where we expect an optimal statistics and signal-to-background 
ratio at space-borne antiparticle detectors p5|, [TOt - We take the experimental sensi- 
tivity to be that of the Pamela experiment after three years of data-taking as our 
benchmark. Finally, the average differential antideuteron flux has been computed in 
the 0.1 < < 0.25 GeV range, where Tjy stands for the antideuteron kinetic energy 
per nucleon, and compared to the estimated GAPS sensitivity for an ultra-long du- 
ration balloon-borne experiment [^] (see Ref. |71| for an updated discussion of the 
role of antideuteron searches in DM indirect detection) . 

In Fig. |, we show various direct and indirect DM detection rates for tuq = mi/2 = 300 
GeV, with Aq = 0, tan (3 = 10 and fx > 0, while M3 is allowed to vary. The M3 value 
corresponding to the mSUGRA model is denoted by a dashed vertical line, while the 
LM3DM scenarios for < and > with ^^_^h'^ = 0.11 are denoted by dot-dashed and 
dot-dot-dashed vertical lines, respectively. The dotted lines correspond to the sensitivity 
level of each of these experiments; i.e., the signal is observable only when the model 
prediction is higher than the corresponding dotted line. While the minimum sensitivity for 
the direct detection rates in frames b) - f) refers to the minimum magnitude of the signal 
that is detectable (and hence independent of the LSP mass), the smallest detectable cross 
section shown by the dotted curves in frame a) depends on the value of m^^ . 

In frame a), we plot the spin-independent neutralino-proton scattering cross section. 
We see that as M3 is decreased from its mSUGRA value, the neutralino-proton scattering 
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mg=300 GeV, m^/2=300 GeV, tan|3=10, Aq=0, ^>0, m=175 GeV 
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Figure 9: Rates for direct and indirect detection of neutralino dark matter vs. M3 for toq = 
mi/2 = 300 GeV, with tan/? = 10, Aq = 0, ^ > 0. Frames c) -f) show the ratio of indirect detection 
rates compared to the mSUGRA modeL In this figure, we adopt the adiabaticaUy contracted N03 
distribution for halo dark matter for our projections of the reach of the various experiments. 



cross section rises almost two orders of magnitude to a value above 3 x 10^ pb, which 
should be detectable by CDMS2, and certainly at stage 3 detectors. A similarly large rate 
is attained for < 0, as shown in the left-hand side of frame a). This frame merely reflects 
the well-known result that MHDM has rather large neutralino-proton scattering rates, as 
is typified by the HB/FP region of the mSUGRA model. The value of asi{Zip) is further 
enhanced by the lowered squark masses of the LM3DM senario. 

In frame h) , we show the flux of muons from neutralino pair annihilations in the core of 
the Sun. The expected muon flux is below the reach of IceCube in the mSUGRA framework, 
but increases over two orders of magnitude into the observable range for IceCube as a result 
of the increased higgsino content of the LSP and the reduced squark mass when the relic 
density is in agreement with the WMAP measurement as in the LM3DM model. 

In frames c), rf), e) and f) we show the flux of photons, positrons, antiprotons and 
antideuterons, respectively. The results here are plotted as ratios of fluxes normalized to 
the corresponding mSUGRA point, in order to give results that are approximately halo- 
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m^/2=''005 GeV, M^<m^^^, tanp=10, A^=0, |x>0, m=175 GeV 
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Figure 10: Rates for direct and indirect detection of neutralino dark matter vs. mo for the 
mi/2 = 1005 GeV slice of the plane in Fig. ^& together with the projected sensitivities of the various 
experimments, assuming the adiabatically contracted N03 distribution for halo dark matter. 



model independent. Also shown by the horizontal lines are the expected experimental 
reaches, as obtained by using the Adiabatically Contracted N03 Halo model|42|. The rates 
for indirect detection via observation of halo annihilation remnants are typically low for 
bino-like DM as in the mSUGRA model. However, when |r3| is reduced until the measured 
CDM relic abundance is achieved, these halo annihilation signals all jump by factors of 
100-200, and are much more likely to be observed by various gamma ray and antimatter 
detection experiments. We should, however, keep in mind that this conclusion is sensitive 
to our assumption of the DM halo profile. 

The detectability of the same signals, but this time for the slice of the mi/2 = 1005 GeV 
slice of LM3DM parameter space we considered in Fig. |^ is illustrated in Fig. |l^. The cross 
sections or the expected fluxes are absolutely normalized, rather than to any particular 
mSUGRA model. The most striking feature of this figure is the rather sudden increase 
(around mo ~ 500 GeV) of the signal as mo increases from low values where the LSP is 
bino-like to high values where the required low value of leads to a significant higgsino 
component in the LSP. 
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We present an overview of the reach of direct and indirect dark matter search tech- 



niques in Fig. 11 and Fig. 12. In both figures, we show the boundary of the region where 



there will be an observable signal in the corresponding dark matter detection channel: the 



signal will be detectable in the direction indicated by the arrows. In Fig. 11 we show re- 
sults for those detection techniques for which the reach is mainly determined by the local 
dark matter density and average circular velocity at the Sun location (here assumed to 
respectively be pDM{ro) = 0.38 GeV cm~^ and U(ro) = 221 km s~^) and is relatively less 
sensitive to the dark matter halo profile, namely direct detection and the fiux of energetic 



neutrinos from the center of the Sun originated by neutralino annihilations. In Fig. |12 
we collect, instead, those quantities whose dependence on the details of the dark matter 
halo is more critical, namely antimatter fiuxes from neutralino annihilations in the galactic 



halo. For the latter, we adopt the Adiabatically Contracted N03 Halo model|42|. We do 
not show the sensitivity contours for GLAST, since with this choice for the dark matter 
halo profile all of the parameter space is within the reach of the space-borne gamma-ray 
detector. 

We notice that while the prospects for Stage-2 detectors do not look particularly 
promising in this context, with the possible exception of small regions at low neutralino 
masses and with light scalars, next generation Stage- 3 detectors will have a sensitivity 
which we estimate to be able to cover most of the parameter space of the models under 
consideration. The sensitivity of IceCube will be instead critically dependent on the hig- 
gsino fraction of the lightest neutralino, which not only controls the pair annihilation cross 
section, but, more critically, sets the spin-dependent neutralino- nucleon scattering cross 
section. The reach contour we find follows in fact quite closely the LSP higgsino fraction. 
However, close to the boundary of the LEP excluded region at the bottom of frames a) and 
b), the neutralino mass is so low that the annihilation neutrinos are not energetic enough 
to give the required muon fiux above the IceCube detection threshold. 

Turning to antimatter searches, we assess here the Pamela sensitivity for primary 
antiproton and positron fiuxes following the approach of Ref. fi^, evaluating the projected 
total aiid demanding a statistically significant (at 95% C.L.) excess over the estimated 
background.^ In Fig. 12, we shade in grey those regions which are already statistically 



excluded at the 95% C.L. by current data on the antiproton fiuxes, when combined with 
an independently estimated secondary and tertiary background [|^. We stress though 
that this exclusion is very sensitive to the assumed halo profile, so is by no means a 
rigorous bound. In general, we find the most promising antimatter search technique will 
be the antiproton channel; remarkably enough, within this scenario we expect a signal at 
space-based antiprotons searches for neutralino masses as large as 0.5 TeV. A low-energy 
antideuteron signal is also expected at GAPS even on a balloon-borne experiment in a 
quite large portion of the parameter space. Finally, a significant positron signal can be 



^Specifically, we evaluate the Pamela sensitivity on the {mo, 7711^2) planes using the approach of Ref. 
psl , where the authors evaluate a prospective taking into account a background independently computed 
with the Galprop package and an estimated energy binning j?^]. This approach is more accurate than 
looking at a single energy bin, although it has been checked that the two approaches are in reasonable 
agreement. 
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Figure 11: Projections for the reaeli of Stage-2 and Stage-3 direct detection experiments, and of 
IceCube in tlie {mo,mi/2) planes of tlic LM3DM model witli r3 < (left) and > (right). We 
shade in red regions where the stau is lighter than the LSP, where electromagnetic gauge invariance 
is spontaneously broken, or where the lightest chargino mass is not compatible with the LEP-II 
bound. The arrows denote the regions where the signal should be detectable. 



only marginally reconciled with current constraints from antiprotons, and might take place 
in the low scalar masses portion of the planes under scrutiny here. 

We notice that in general the anti-particle sensitivity contours we obtain trace the 
higgsino fraction of the LSP: for very small values of mo, the sensitivity drops at lower neu- 
tralino masses even though the antimatter primary fluxes approximately scale as (av) /rn- , 
because the LSP becomes more bino-like. Notice also the corridor around mi/2 = 400- 
450 GeV where the current anti-proton searches do not exclude the model. For smaller 
values of mi/2) the LSP mass is small and the anti-proton flux is too large. This flux falls 
below the 95% CL limit until the LSP mass becomes large enough so that annihilation 
to tt becomes kinematically allowed, once again yielding a large anti-proton rate. Finally, 
when mi/2 gets larger, the LSP mass increases, and the anti-particle rate once again drops 
below the current experimental limit. 

We stress that the exclusion limits shown in frames c) and d) as well as the projections 
for the reach via anti-particle and gamma ray searches are sensitive to our assumption of 
the adiabatically contracted N03 halo profile which yields considerably higher values for anti- 
particle, and especially, gamma ray fluxes. Assuming different but equally viable galactic 
DM halo distributions |4^, |l^ modifies this conclusion. Until the halo profile can 

be independently determined, we believe that exploration of independent signals even in 
these "halo-profile-dependent excluded regions" should continue. 
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Figure 12: The same as Fig. m, but for antimatter searclies at PAMELA and GAPS. GLAST 
will have the sensitivity to probe the entire plane. These projections are sensitive to our choice of 
the adiabatically contracted N03 halo model for the distribution of the galactic dark matter. The 
regions shaded in grey are excluded by current antiproton data, again assuming the adiabatically 
contracted N03 halo model, but not necessarily for other equally viable halo profiles. 



4. LM3DM at colliders 
4.1 Fermilab Tevatron 

In SUSY models with gaugino mass unification, the ratio of weak scale gaugino masses is 
typically found to be Mi : M2 : M3 ~ 1 : 2 : 7. In the mSUGRA model, since ~ M2 
and rrig ~ M3, the bound on chargino masses > 103.5 GeV from LEP2 implies as 

well that nig ~ 350 GeV. Since a 400 GeV gluino is typically beyond the reach of Fermilab 
Tevatron experiments, this leaves a relatively tiny window for a gluino discovery at the 
Tevatron, at least within the mSUGRA framework[74|.^ 



However, in the case of LM3DM, the gluino mass can be much reduced relative to the 



value of rn^^^. The situation is illustrated in Fig. 13, where we plot contours of rrig in 



the same mo vs. mi/2 plane as in Fig. |5| where has been dialed to low enough values 
that the value of ^^_^h'^ = 0.11 everywhere. Of importance here is that the gluino mass 
immediately adjacent to the blue-shaded LEP2 excluded region reaches values below 200 
GeV, which is surely accessible to Fermilab Tevatron searches. In fact, for niQ = 1500 GeV, 
mi/2 = 250 GeV, we find the LM3DM spectrum labelled as LM3DM3 in Table ^ where 
rrig = 183 GeV is consistent with LEP lower limits on the chargino mass. In this scenario, 
the cross section for pp — > gg at the Fermilab Tevatron with ^/s = 2 TeV is ~ 21.5 pb, 
so that for 1 fb^^ of integrated luminosity, we already expect in excess of 20K gg events! 



®For m-g = 400 GeV and = 2mg, we find a{pp gg) = 27.8 fb at ^ = 2 TeV. 
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Figure 13: Contours of nig in tlie rno vs. plane for tan/3 = 10, = 0, /i > and a) 

mSUGRA model, b) A/3 < LM3DM and c) M3 > LM3DM. 



Thus, considerable portions of the mo vs. mi/2 parameter plane that could not be probed at 
LEP in the LM3DM scenario should be accessible to present day Tevatron SUSY searches! 

Once the gluino pairs are produced, it is important to examine their decay modes. 
Inspection of the Isajet decay table for point LM3DM3 shows the surprising result that 
in this region of parameter space, gluino decays via three body modes g — > qqZi and 
9 ~^ QQ'Wj are suppressed, and that in fact the loop decays g — > Zig are dominant! This 
large enhancement of the radiative decay relative to three body decays is not hard to 
understand. It has long been known |75t| that, because of the tracelessness of the diagonal 
generators of the electroweak gauge group, for degenerate squarks, all contributions from 
the gaugino components of the neutralino cancel in the amplitude for the decay g gZi. 
As a result, unless tan/3 is very large, top squark loops completely dominate the radiative 
decay amplitude. The usual three-body decays of gluinos, on the other hand, receive 
significant contributions from the gaugino components of neutralinos, and indeed dominate 
these decays unless tan P is very large. For very heavy squarks, neglecting bottom quark 
Yukawa couplings and phase space effects, we find that 

r(g ^ gZ^) _ 12 
T{g ^ qqZi) ^ 

where, in the notation of Ref . p^] , ^ ~ and ~ are couplings of the i*^ neutralino (gaugino 

components) to the quark-squark system, and v'f^ is the component of the higgsino field 
(that couples to up type (s)fermions) in this neutralino. In deriving ([4.1|), we have assumed 
that the v^^'^ <C max{v^r^^'^ ,v'f^'^), i.e. that the neutralino is mainly gaugino-like. In many 
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LM3DM: -M^ < m^,^ , tan^=10, =0, \i >0, m, =175 GeV 



LM3DM: < m^,^ , tanp=10, 4„ =0, n >0, m, =175 GeV 
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Figure 14: Contours of BF{g Zig) (summed over i = 1 — 4) in tlie mo vs. miji plane for 
tan /5 10, Ao 0, > and a) M3 < LM3DM and b) M3 > LM3DM. 



models, where light neutralinos have only small higgsino components, gluino radiative 
decays have very small branching ratios because of the factor in This same 

factor is, however, precisely the reason for the large branching fraction in the LM3DM3 
scenario that we have been discussing. 

To bring home this point, we show in Fig. 14 the g — > Z^g branching fraction contours 
summed over i = 1 — 4 in the mo vs. mi/2 plane of Fig. ^. We see in the lower right region 
where gluino masses are quite light that the cumulative gluino loop decay branching exceeds 
Thus, in this region, gg production events with g 



Z\g will give rise typically 
to dijet+ii^™'^^ events, much like squark pair production when rriq < nig. However, the 
g Zig decays are not only into gZi , but also have large rates into Z2g and Z2,g, as shown 
in Fig. 15, where these branching fractions are plotted versus mo for fixed mi/2 = 300 



GeV and other parameters as in Fig. |T^. 

The Z2 and Z^, which are produced either directly or via gluino cascade decays will 
likely decay via three-body modes which, if rriq is large enough, are dominated by Z ex- 
change. The m^^ — m^^ mass gap is shown in Fig. Since |^| is typically quite low, 
and the lighter Zi are mixed higgsino states, this mass gap varies in the 30-70 GeV range 
when Zi is MHDM, so that Z2 UZi (and also frequently Z3 UZi) occur all over the 
MHDM portion of the LM3DM parameter space. In this case, one or possibly even two 
distinct m{£^£~) mass edges should be apparent if enough sparticle pair production events 
are generated. The m{i^i~) mass edges of course are renown for being the starting point 
for sparticle mass reconstruction at hadron colliders f?^ . 

4.2 CERN LHC 

When comparing the LM3DM scenario to the mSUGRA model, we have found that for 
a given point in the mo vs. mi/2 plane, gluino and squark masses are quite suppressed 
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Figure 15: A plot of BF{g Z^g) (for i = 1 - 4) versus toq for mi/2 = 300 GeV, Aq ^ 
tan/3 = 10, /i > where at each point M3 has been dialed so that = 0.11 



relative to slepton masses and somewhat suppressed relative to chargino masses. The 
upshot is that sparticle production at the LHC should be even more dominated by gluino 
and squark production (compared to other sparticle production reactions) than in the case 
of mSUGRA. Indeed, for the points listed in Table |l] we find (j{pp ggX) = 2, 84, 31 and 
4760 pb , for mSUGRA, LM3DM1, LM3DM2 and LM3DM3, respectively. The signatures 
from gluino and squark pair production at the LHC from LM3DM will consist of various 
multi-jet plus multi-lepton plus i?™^^ events as in mSUGRA [ 77 1. However, in the lower 
right portion of the mo vs. plane, the SUSY events should consist mainly of dijet 

_l_^miss gvents when g Zig, with additional jets and opposite sign/same flavor (OS/SF) 
dilepton events coming from g — > Z2g or Z^g followed by Z2 or Z3 decay. In particular, 
same sign dileptons, which are somewhat characteristic of gluino pair production |7^, will 
be relatively suppressed when the gluino loop decays are dominant. 

The LHC reach for SUSY in the mSUGRA model has been calculated in Ref. (7|, |0| 
assuming 100 fb~^ of integrated luminosity. The ultimate reach of the LHC is mainly 
dependent on the value of rrig and niq, and not so dependent on their particular decay 
modes. Thus, we may translate the mSUGRA reach results into contours into the rriq vs. nig 
plane, and then convert this contour into a reach contour in the mo vs.mi/2 plane of the 



LM3DM scenario. The translated reach contour is shown in Fig. [17| for the same parameters 
as in Fig. |5|6) . The reach contour on the left portion of the plot in the BDM region mainly 
follows the rUg ~ niq ~ 3 TeV contour. The upper left increase in reach is due to the 
sliver of MHDM region at low mo and large mi/2 shown in Fig. where is relatively 



- 22 - 



mSUGRA: tanp=10, A,, =0, \i >0, m, =178 GeV 
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Figure 16: Contours of m^^ — m^^ mass gap in tlie mo v.s. rni/2 plane for tan/3 = 10, = 0; 
/i > in the case of a) the mSUGRA model, 6) M3 < LM3DM and c) M3 > LM3DM. 

reduced. At intermediate niQ, neutralino annihilation is assisted by the A resonance, and 
so higher values are found to match the WMAP constraint (recall in this region the Zi 
is bino-like). On the right hand side of the plot for mo > 1-2 TeV, we are in the MHDM 
region with lower values, and consequently lighter gluino and squark masses. Thus, the 
reach is increased. At very large mg values, rriq > rrig, and here the LHC reach extends 
only out to nig ~ 2.7 TeV. We should mention that since light neutralinos and chargino 
have significant higgsino component in the LM3DM scenario, gluino decays to the third 
generation quarks will be enhanced exactly as in the FP/HB region of the mSUGRA model, 
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so that 6-jet tagging may improve the LHC reach by ~ 10 — 15% beyond that shown in the 

LM3DM: < m^^^ , tanp=10, Ag =0, \i >0, m,=175 GeV 



figure. ||81 1 




1.8 2 
m„(TeV) 

Figure 17: Reach contours for the CERN LHC with 100 fb^^ of integrated luminosity and for a 
y/s = 0.5 and 1 TeV linear e~^e~ collider in the mo vs. mi/2 plane for tan/3 = 10, Aq — 0, > 
and where has been reduced such that ^^_h^ =0.11 at every point in the plane. 

For SUSY searches at the CERN LHC, Hinchhffe et al. have pointed out|Q] that an 
approximate value of rriq or rrig can be gained by extracting the maximum in the M^ff 
distribution, where Me// = + Erijet 1) + Erijet 2) + Erijet 3) + Erijet A). Their 
analysis will carry over to much of the LM3DM scenario,"^ as well as in models with gaugino 
mass unification, so that the approximate mass scale of strongly interacting sparticles will 
be known soon after a supersymmetry signal has been established. 

In mSUGRA, a dilepton mass edge should be visible in SUSY signal events only if 
mi/2 ~ 250 GeV (where Z2 Z\it is allowed) or if two body Z2 U, U decays are 
allowed. In the case of LM3DM, as with MWDM[||] and BWCA DM[|3|, the dilepton 
mass edge should be visible over almost all parameter space. We illustrate the situation 
for three of the case studies listed in Table |l[^ The first case, labeled mSUGRA, has 

^An exception is the low mi/2, large mo region of Fig. |^ where dijet production is instead dominant. 

*In this study, a toy detector simulation is employed with calorimeter cell size A77 x A<jf> = 0.05 x 0.05 
and —5 < < 5. The hadronic energy resolution is taken to be 80%/-\/S for \ti\ < 2.6 and 100%/\/S 
for \ri\ > 2.6. The electromagnetic energy resolution is assumed to be 3Vo/y/E. We use a UAl-like jet 
finding algorithm with jet cone size 7? = 0.5 and py' > 25 GeV. We also require that \rii\ < 2.5 and 
|77j| < 3. Leptons (es or ^s) have to also satisfy pl^^'°" > 10 GeV. Leptons are considered isolated if the 
visible activity within the cone AR < 0.3 is Ei55^'^"° < 2 GeV. The strict isolation criterion helps reduce 
multi-lepton background from heavy quark (especially tt) production. 
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mg = rni/2 = 300 GeV, with Aq = 0, tan/3 = 10 and fi > 0. In this case, gg, gq and qq 
production occurs with a combined cross section of about 12 pb, while the total SUSY cross 
section is around 13.4 pb (the additional 1.4 pb comes mainly from -ino pair production and 
-ino-squark or -ino-gluino associated production). The case of LM3DM1, with M3 = 150 
GeV, has a total cross section of 215 pb. The case of LM3DM2, with slightly heavier 
squark and gluino masses, has a total production cross section of 101 pb. The special case 
of LM3DM3, which should be accessible to Tevatron searches via its light gluino, has a 
total SUSY cross section at LHC of 3744 pb. 

We have generated 750K LHC SUSY events for the cases LM3DM1 and LM3DM2 using 
Isajet 7.73, and passed them through a toy detector simulation as described above. We 
adopt cuts which are similar to those of LHC point 5 of the study of Hinchliffe et a/. p2|| , 
which efficiently select the SUSY signal while essentially eliminating SM backgrounds: 
^miss y rnax{WO GeV, 0.2Me//), at least four jets with Et > 50 GeV, where the hardest 
jet has Et > 100 GeV, transverse sphericity St > 0.2 and Meff > 800 GeV. 

In these events, we require at least two isolated leptons, and then plot the invariant 



mass of all same flavor/opposite sign dileptons. The results are shown in Fig. 18. In the 
case of the mSUGRA model, frame a), there is a sharp peak at m{i~^i^) ~ Mz, which 
comes from Z2 — > ZiZ^ decays where Z2 is produced in the gluino and squark cascade 
decays. In the case of LM3DM1 in frame b), we clearly see a continuum distribution with 
a mass edge at m{i^£^) < m^^ — m^^ = 58.5 GeV. We also see events beyond this edge 
along with a peak at Mz- In this case, m^^ — m^^ =90.3 GeV is within Tz of Mz and we 
would expect that dileptons from Z3 Z\ll decays have their mass sharply peaked just 
below Mz- This peak would also be populated by Z bosons produced via Z4 ZiZ^ or 
W2 — > WiZ^ decays. The cross section plotted here is ~ 188 fb, which would correspond 
to 19K events in 100 fb~^ of integrated luminosity (the sample shown in the figure contains 
just 646 events). In frame c)- with a cross section of ~ 207 fb (but just 1550 actual entries)- 
once again we see the Z^ peak from decays of heavier charginos and neutralinos to the Z 
boson, together with a mass edge at m{i^i^) < m^^ —m^^ = 61 GeV, and a continuum in 
between, presumably mainly from chargino pairs in SUSY events. In both these LM3DM 
cases, the m^^ ~ ''^Zi edge should be very precisely measurable. It should also be 

clear that this edge is inconsistent with models based on gaugino mass unification, in that 
the projected ratios Mi : M2 : M3 will not be in the order 1 :~ 2 :~ 7 as in mSUGRA. 
Although the Z2 — Z\ mass edge will be directly measurable, the absolute neutralino and 
chargino masses will, as usual, be more difficult to extract at the LHC. 

4.3 Linear e^e~ collider 

The reach of the CERN LHC for supersymmetric matter is determined mainly by rriq 
and m^, which depend on mo and M3. In contrast, the reach of the ILC for SUSY is 
largely determined by whether or not the reactions e~ — > M^^^'W^]" or e^e~ — > i^'t~ 



are kinematically accessible [83]. For instance, chargino pair production is expected to be 
visible if -^s > 2m^^. The value of J^^y^ depends mainly on M2 and ^i- Thus, in the 
LM3DM scenario where M\ and M2 take values similar to mSUGRA, but where /i is quite 
small, the reach of the ILC in the mg vs. mi/2 plane via chargino pair production will be 
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Figure 18: Distribution of same flavor/opposite sign dileptons from SUSY events at the CERN 
LHC from a) mSUGRA, b) LM3DM1 and c) LM3DM2 cases as in Table ||. 



enhanced relative to the case of the mSUGRA model. Since slepton masses are relatively 
unaffected by lowering M3, the ILC reach for slepton pair production will be similar to 
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LIVI3DM: < m,^^ ' tan^=10, Ag =0, h >0, m, =175 GeV 




Figure 19: Reach of the ILC for SUSY in the LM3DM scenario where AI3 is lowered until 
57^^/1^ — 0.11 at every point in the mo vs. mi/2 plane. We show the ILC reach for ^/s = 0.5 TeV 
and 1 TeV via the kinematic limits for W^W^ , f+f~ and iiii production. 



the mSUGRA case. In addition, squark masses are relatively suppressed in the LM3DM 
scenario, especially the top squark, so that there will be a non-trivial reach of the ILC for 



titi production. The situation is illustrated in Fig. 19 where we show the ultimate reach of 
the ILC in the tuq vs. mi^2 plane for tan /3 = 10, j4o = 0, /x > and rrit = 175 GeV. We have 
dialed M3 at every point to give ^^^^^^ = 0.11, in accord with the WMAP observation. The 
reach of ILC with ^/s = 0.5 TeV is denoted by dashed contours, and extends to mi/2 ~ 500 
GeV, while the corresponding reach within the mSUGRA framework with gaugino mass 



unification extends to ~ 320 GeV[83|. The reach of ILC with ^/s = 1 TeV extends to 
mi/2 ~ 1 — 1.2 TeV, compared with the mSUGRA value of mi/2 ~ 600 GeV. The combined 
reach (from chargino and selectron production) of the ^/s = 0.5 and 1 TeV ILC relative 
to the LHC are shown in Fig. |l^. The LHC reach is always larger than that of the ILC, 
primarily because of the relative reduction of gluino and squark masses in the LM3DM 
framework. 

The distinguishing feature of LM3DM is that the small fi parameter gives rise to a 
rather light spectrum of the two charginos and all four neutralinos. Thus, many more -ino 
pair production reactions are likely to be accessible to a linear e^e~ collider than would 
occur in the mSUGRA model. As an example, we show in Table ^ the various SUSY cross 
sections in fb which are accessible to a = 0.5 TeV machine for the four case studies listed 
in Table ||. While only WiWi, Z1Z2 and Z2Z2 production are possible in the mSUGRA 
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Table 2: Cross sections in fb for e+e SUSY processes at a y/s = 0.5 TeV linear collider, for 
the four case studies listed in Table 1. 



model, for the LM3DM scenarios all ten reactions listed are accessible at a ^/s = 0.5 TeV 
linear collider, although some of these have very low rates. It does appear though that 
every chargino and neutralino is produced via some reaction with cross section exceeding 
10 fb. Detailed studies of the chargino-neutralino sector along the lines discussed in Ref. 
\8§ should be feasible within the LM3DM scenario. 

Another feature of LM3DM relevant to linear e^e~ colliders is that the relatively low 
squark masses which are expected in this scenario means that squark pair production is 
more likely to be possible, especially for a ^/s > 1 TeV machine. In most cases, the decay 
q —> gq \s kinematically allowed, so that gluino production might be studied in the e+e" 
environment, since they can be produced via the cascade decays of the heavier squarks. In 
this case, very precise determination of squark and gluino masses may be possible if the 
end point of the energy spectrum of the primary quark jet in the decay q ^ qg can be 
identified. 



5. Summary and conclusions 

If we identify the relic density of CDM ( |l.lD , with that of thermally produced LSPs of an 
R-parity conserving SUSY model, the WMAP measurement serves as a stringent constraint 
on any SUSY model. It is then interesting to explore the ramifications of this measurement, 
both for collider searches for supersymmetry, as well as for direct and indirect searches for 
DM at non- accelerator facilities. It is also necessary to explore just how robust these rami- 
fications are to changes in the underlying SUSY framework that do not alter the successful 
prediction for the LSP relic density. In previous studies, we explored how the WMAP 
CMD constraint could be satisfied if (1) we relax the (phenomenologically unnecessary) 
assumption that the Higgs scalar mass parameters unify with sfermion mass parameters at 
high scales, and (2) if we relax the assumption of the unification of gaugino masses, and 
allow the ratio M1/M2 (which controls the composition of the LSP) to float freely. In this 
paper, we study the implications of what we dub as the low M3 DM (LM3DM) model, 
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which is essentiahy the paradigm mSUGRA framework, except that the SU{3) gaugino 
mass is allowed to adopt any value. Following earlier studies |25, 26| we find that, for 
mo values not hierarchically larger than mi/2, the value of |/i| is reduced when the GUT 
scale gluino mass parameter IM3I < mi/2, and MHDM or BDM solutions become viable 
for essentially all values of model parameters. 

The sizeable higgsino component of MHDM implies enhanced detection rates in on- 
going, planned and proposed experiments searching for DM relative to the bino LSP case 
more typical in mSUGRA; see Fig. 11 and Fig. 12. While direct searches at stage 2 detectors 
such as CDMS2 can explore only a relatively limited portion of the parameter space, the 
entire parameter space should be explorable at the proposed stage 3 detectors typified by 
the SuperCDMS or 1-ton XENON experiments. Indirect searches via the detection of hard 
muon neutrinos from the core of the sun should also be possible at IceCube over much of 
the model parameter space. Experiments looking for anti-particles and gamma rays from 
the annihilation of neutralinos in our galactic halo should also be able to detect signals 
from MHDM. These projections should be interpreted with care because they are sensitive 
to the precise distribution of the DM in the galactic halo. 

By comparing detection rates in direct and indirect search experiments, it is possible 
to qualitatively distinguish the MHDM scenario from scenarios where the dark matter is 
bino-like as in mSUGRA (either in the bulk region or in the Higgs resonance region) or 
with bino-wino coannihilation yielding the WMAP value, or mixed with the wino p3|| . 
Experiments at colliders will be able to provide additional evidence in favor of one or the 
other of these scenarios. 

MHDM, on the other hand, occurs in a variety of models. It may occur for very large 
values of mo in the HB/FP region of the mSUGRA model, in non-universal Higgs mass 
(NUHM) models where the GUT scale Higgs mass parameters are equal but larger than 
the corresponding sfermion parameters or in more general NUHM models, or, as we have 
just seen, in the LM3DM models. Distinction between these various MHDM scenarios is 
only possible via examination of the properties of other sparticles which are accessible via 
collider searches for SUSY. 

The main distinguishing feature of the LM3DM model is that the ratio of coloured 
sparticle masses to those of colour singlet sleptons, charginos and neutralinos is smaller 
in the LM3DM model than in most other models. This clearly favours SUSY searches 
at hadron colliders such as the Fermilab Tevatron or the CERN LHC vis a vis searches 
at electron-positron colliders. For instance, while the LEP lower limit on the chargino 
mass greatly restricts the potential of the Fermilab Tevatron to discover gluinos within the 
mSUGRA framework (or, for that matter, in any framework with unification of gaugino 
masses), Tevatron searches for gluinos may yet lead to the discovery of SUSY if SUSY is 
realized as in the LM3DM model\85\. In this case, experiments at the LHC will have a 
reach much larger than that of even a TeV linear collider. Despite this, experiments at 
the linear collider will play a big role in elucidating the physics and allowing us to zero in 
on the underlying scenario. Since is comparable to the weak scale electroweak gaugino 
masses, it is likely that all charginos and neutralinos will be accessible and their properties 
measured at a TeV linear collider. In this case, it will be possible to directly determine Mi, 
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M2, jJL and tan [5. Combining these with the determination of nig that should be possible at 
the LHC, we should be able to determine that the GUT scale gluino mass is smaller than 
the corresponding electroweak gaugino masses.^ If we are lucky, the top squark and perhaps 
even other squarks, may be kinematically accessible. In this case, the gluino may also be 
accessible as a decay product of the squarks, and true bottom-up sparticle spectroscopy 
would be possible. 
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